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I. INTRODUCTION 



Bianchi cosmological models have long been used as toy models in both classical general 
relativity and quantum gravity. The small but non-trivial number of degrees of freedom they 
possess makes their dynamics interesting from the point of view of shedding light on features 
of the full theory, even though it reduces to that of a relativistic particle in a low-dimensional 
space. These models are homogeneous cosmologies, in which the symmetries are imposed 
on the canonical action or Hamiltonian before the dynamics is obtained. This procedure 
is, in general, not equivalent to imposing the symmetries on the full dynamical equations, 
but in a classical theory this can only potentially give some extra solutions with respect to 
the symmetric ones for the full theory, and is not a problem 0]. Thus, the classical Bianchi 
models can be very fruitfully used in exploring issues such as the evolution of anisotropics in 
simplified versions of the early universe, or the chaotic nature of the evolution of the metric 
close to the initial singularity. 

In the quantum theory, the difference between any physical prediction of a minisuper- 
pace model and the "full theory" may be more substantial f5], but it is generally believed 
that some qualitative features of the models are still useful indications of what to expect 
from the full theory when different approaches to its quantization are used. In particular, 
minisuperspace models are time-reparametrization invariant, thus exhibiting the much cele- 
brated "problem of time" they represent then a fruitful playground in which to explore 
quantization methods, in the hope of getting some insight on how to approach this problem 
in full quantum gravity. It is not surprising then that minisuperp^aces have been analyzed 
from the Euclidean path integral approach to quantum gravity [4], consistent histories |^ 
and refined algebraic quantization [6] perspectives. 

It is in this spirit that we study here coherent states of minisuperspace models. The 
aim is to understand the role of coherent states in defining semiclassical states for time- 
reparametrization invariant systems. In recent years, an increasing amount of work has been 
devoted to understanding the phenomenology of semiclassical states in quantum gravity, by 
which here we mean that sector of the physical states in which the gravitational degrees 
of freedom are sharply peaked around values satisfying the classical constraint equations, 
because it is expected that these states will possess quasi-classical properties in certain 
regimes 0. The study of the semiclassical limit of quantum cosmological models is of 
course not new; however, most of the treatments that we are aware of deal with the WKB 
method, where the wave function is assumed to be of a special form and some component 
of it is shown to satisfy the Hamilton- Jacobi equation (for a review, see Ref jjSj). 

In this paper, for simplicity, we consider vacuum Bianchi type I models, whose classical 
description can be reduced to that of a relativistic particle moving in a 3-dimensional, linear 
configuration space, subject to the Hamiltonian constraint of general relativity, which in this 
case reduces to the free particle Hamiltonian. These models have been extensively studied, 
but the type of questions we are concerned with here, however, have been only recently 
considered within the context of constrained systems in Ref 0], where a general framework 
for the construction of semiclassical states for a class of constrained theories is discussed. 
Here, we follow much of that approach in spirit, but the constraint and physical observables 
that arise in these Bianchi models are not of the same type and we have to adapt the general 
constructions to the present situation, and the fact that we are dealing with cosmological 
models allows us to ask additional physically interesting questions. 

Despite the simplicity of the classical model, the quantum counterpart is non-trivial, be- 
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cause of the presence of the constraint. First, we consider coherent states in the kinematical 
Hilbert space, peaked around a point on the constraint surface in phase space but without 
taking the quantum constraint into account, and characterized by a set of "squeezing" pa- 
rameters. These states are only meant to be approximations to physical coherent states, 
i.e., solutions to the quantum constraint, which we find by applying the group averaging 
procedure of the Refined Algebraic Quantization program (RAQ) p to the kinematical 
coherent states; one expects these to be the ones that capture the "essential" features of 
semiclassical Bianchi I models. Finally, we consider the reduced phase space, described 
by true Dirac observables, and define coherent states peaked around reduced phase space 
points; one might call these "intrinsic" coherent states. The main purpose of this article 
is to compare these three classes of states, in terms of describing semiclassical situations, 
assess how good the kinematical coherent states are as approximations to the physical ones, 
and use the "intrinsic" ones to get an insight into evolution issues. 

The structure of the paper is the following. In Sec. Owe recall some basic notions about 
Bianchi models. In Sec. IIIII we recall the notion of coherent states, and discuss what we 
mean by semiclassical states for a system. Sec. II VI considers the kinematical coherent states 
for Bianchi I models, and in Sec. |3the Physical coherent states. In Sec. IVII we study the 
time evolution of the "intrinsic coherent states," and in particular some estimates of the 
time it would take for a coherent state to lose its peakedness properties. We end with a 
discussion in Sec. IVIII Throughout the paper we set c = G = 1, but keep h dimensionful. 

II. PRELIMINARIES ON BIANCHI MODELS AND THEIR OBSERVABLES 

In this section, we will provide some background information on Bianchi models, the vari- 
ables used to describe them, and their physical observables; readers who are familiar with 
these models should be able to skim through it quickly. 

Bianchi models are spatially homogeneous cosmologies: they describe spacetimes which 
can be foliated by spacelike hypersurfaces (all diffeomorphic to the "space manifold" S), 
on each of which a group G of isometrics acts transitively, and they are classified by the Lie 
algebra of the group G. 

If one uses three linearly independent left- invariant forms a* on S as a basis for its 
cotangent space, the spacetime metric can be written in the form g = — iV^(t) dt dt + 
qij{t) (T* ® cr-? , where in some Bianchi models the matrix gij can be diagonalized and written, 
in the Misner parametrization as 

qi, = e'"" diag(e2(^++v^^"),e2(^+-v^/3-)^g-4/3+) . (2.I) 

we thus see that -v/det q = e^'^" , and characterizes the volume of E^, while the 

characterize its anisotropy; all three functions take values in (— oo,-|-oo). We will take as 
configuration variables the P\ with i = 0,+,—; the configuration space is thus C = IR^. 
The dynamics of these models then reduces (before imposing the constraints) to that of a 
relativistic particle in 3 dimensions, and is governed by a Hamiltonian constraint of the form 

^ = h^.Pi + (2.2) 

in the vacuum case, where pi is the canonical momentum conjugate to and rjij a flat metric 
on C, with line element r]ij d/S^d/S^ = — (d/3°)^ + (d/3+)^ + (d/3~)^, and the form of the potential 
U depends on the model under consideration. Here we should point out that even when 
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the flat metric resembles that of 3D Minkowski spacetime, it is actually conformally related 
to the De Witt supermetric on C. As with all general relativistic models, the Hamiltonian 
is also a constraint, in that allowed classical phase space points have to satisfy the scalar 
constraint H{pi,f]^) = 0; this is the only constraint for these models, since the metric (|2.ip 
satisfies the vector constraint identically. 

In this paper, we restrict ourselves to begin with to the simplest models, the vacuum 
ones of type I (with homogeneity group U(l)'^, or the additive IR^). In this case, U{P^) 
vanishes identically; this allows us to not only solve the quantum constraint easily, but 
also to conceptually isolate the issue of comparing the different coherent states on a linear 
phase space with a relatively simple constraint, from those related to the structure of phase 
space or the form of the potential. In the canonical quantum theory, following the Dirac 
prescription, states have to be annihilated by the operator version H of the constraint, i.e., 
be invariant under the group e~'^'^ it generates, and it is therefore very advantageous, to use 
the p-representation, because the Hamiltonian is then just a multiplication operator. We 
will also rescale the variables by pi ^ pi/h, which makes the new pi dimensionless (like the 
we can think of this alternatively as a rescaling of the action or Hamiltonian, or as a way 
to give the value of pi in terms of the quantum unit h, which is useful when characterizing 
the extent to which a fluctuation is small or large in the classical sense. For simplicity of 
notation, in the rest of the paper we will use only subscripts, and denote Pi = /3* (that is, 
the index is not lowered with the metric rjij). 

In order to evaluate candidate semiclassical states by comparing fluctuations of observ- 
ables, we need to identify the proper set of physical observables to consider. It is natural 
to look for a sufficient number of Dirac observables, in the sense that their values consti- 
tute a full characterization of a (classical, gauge-invariant) state of the system. Since we 
will be using the p-representation, to simplify the terminology it is convenient to call the 
space C = JR^ of piS our new configuration space. Classically, configuration observables are 
then functions f{pi) on C, and "momentum" observables, specifically linear functions of the 
canonically conjugate Pi, can be thought of as vector fields u onC, since there is a canonical 
assignment of a linear momentum observable P„ everywhere in the phase space F to each 
such vector field; we can then associate to these observables quantum operators / and u 
acting on ipi^p), defined respectively by 

Hp) , :=i(/:„ + idiv^n)^(p) , (2.3) 

where div^u is the divergence of u with respect to some volume element d/x(p) on C. 

We know that the reduced phase space has four dimensions, so we need four independent 
observables {00}^=!^... ^4, in addition to the constraint itself. The first observation is that the 
variables p± are Dirac observables, since the Hamiltonian constraint in the new variables, 

HiPi,p,) = l{-pl+pl+pl), (2.4) 

depends on the pi only. The not so obvious task is that of finding a suitable pair of linear 
momentum observables, or vector fields on C, to represent the conjugate variables, since 
the Pi, for which the CCR with the Pi are of the form [Pi,Pj] = iS^j, are not physical 
observables; in fact, one can easily see that no function of the Pi alone is a constant of 
the motion. As Dirac "momentum" observables for this model we will use instead the 
combinations v± := poP± + p±Po which, together with Vo := p+P- — P-P+ are associated 
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with the vector fields (denoted here by the same symbols) 



d d d d d d . . 

v+=Po^—+P+T^, V-=Pq^+P^ — , Vo=p+- P-^— , (2.5) 

opj^ opo op^ opo op- op+ 



i.e., the two boosts and the spatial rotation that generate the (2+1) Lorentz group 3| on 
the 'minisuperspace' C. 

We may note that in general there is a possible ambiguity in the class of observables, 
consisting in adding to each one a term proportional to the constraint; such an addition 
does not change the value of the observable on the constraint surface. One can implement a 
general strategy to remove the ambiguity but in the present case no such change would 
preserve the property that all observables are linear in the configuration variables. 



III. SEMICLASSICAL AND COHERENT STATES 

In this section, we will specify our criteria for a quantum state to be semiclassical, following 
the general discussion of Ref [3], recall the basic definition of the most common example 
of such a state, a coherent state, and specify the setup and goals for the rest of the paper, 
specializing the discussion to Bianchi I cosmolo gy. (For other references on semiclassical 
states for systems with constraints, see also Refs [13,llJl-) 

Given a point {qi,Pi) in a classical phase space, a complete set of physical observables 
{Oa}a=i,...,4: for the theory (completeness will be the only requirement, the observables do 
not need to generate an algebra; they may be overcomplete), and a pair of parameter values 
(ea, 6a) for each Oa, we will consider a state ip to be semiclassical if for each a 

\{da)ij - Oa{qi,Pi)\ < ea , (A6q)^ < 5^ • (3.1) 

In other words, we are assigning tolerances; the {ea,Sa) give a measure of how close the 
expectation values must be to the classical values, and how narrowly peaked a state must 
be, for the latter to be considered semiclassical. One often takes the Oa to be the canonical 
variables {qi,Pi) themselves, but this is not necessary; in some cases, as we will see, one is 
forced to make other choices, and in general the choice can affect in a crucial way which 
states are considered to be semiclassical. 

If the system has a set of constraints Hj{qi,pi) = (assumed here to be first class), a 
couple of additional points need to be addressed. First of all, the observables need to be 
physical in the Dirac sense of commuting (weakly) with the constraints, {Oa,Hj} ^ 0, but 
they need only be complete in the reduced phase space — if one knows how to characterize 
it. Secondly, if the state does not belong to the physical Hilbert space, we specify the value 
of an additional pair of tolerance parameters (e/, 6j) for each constraint, and require that V/ 

\{Hi)^\<ea, iAHj)^<5i. (3.2) 

An important issue, which in fact underlies the motivation for the current work, arises 
when one has a constrained system for which one prefers, or is forced, to work with kinemat- 
ical states. One may want to know in that case whether the kinematical semiclassical states 
are, in an appropriate sense, good approximations to the physical ones, so that one may 
extract from the former useful information on the latter, for example regarding quantum 
fluctuations. The general formulation of this issue deserves further study; here we will just 
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mention that two possible definitions of wlien T/'kin is a good approximation to ipphy are: (i) 
Tlie two states are semiclassical with respect to the same {qi,Pi), with the same tolerances; 
and the more restrictive definition (ii) For all Oa 



|(0.)phy-(0.)kin| l(AQ^)phy-(AOa)kinl 



|(0,)phy| (A0„) 



2 

phy 



< 1 . (3.3) 



Qualitatively, the first definition simply aims at identifying situations in which the spreads 
of measured values for the Oa in ■j/'km and ipphy are both within the same tolerances, while 
the second one may be suited for situations in which one would like all quantities which can 
be calculated from the spread of values of Oa with the two states to nearly coincide. 

It would be natural to ask at this point whether the values of the tolerances {ea, Sa) can 
be arbitrarily chosen. Although the definition above is consistent for arbitrary (positive) 
values, if the 6a are too small no state will satisfy the conditions. For example, in the case 
of a pair of conjugate variables (p, q), choosing 6p and 6g such that 6pSg < fi/2 would be too 
restrictive; in general, quantum uncertainty relations will translate into some (macroscopi- 
cally) mild conditions that set lower bounds for combinations of (5a 's. The in principle can 
be arbitrary, but in practice they may be restricted by the way the states are constructed, 
as mentioned below, and by their interpretation. A physically motivated choice would be to 
set the value of each ea to be close to, or smaller than, that of the corresponding 5a, with 
the latter being smaller than the experimental resolution for Oa- 

We should also comment here that it may seem like the introduction of the parameters 
ea is unnecessary, since one can usually choose the states in such a way that the expectation 
values of the Oa agree with their classical values. This is correct, in principle; however, 
some of the physical states we will use are constructed (using group averaging) out of states 
with the right expectation values for our choice of physical observables, but do not have this 
property themselves. Besides, it may sometimes simplify calculations to choose states which 
don't have exactly the classical expectation values for the chosen Oa- 

Coherent states are among the most commonly used states with a semiclassical interpre- 
tation; they are relatively easy to work with, are often associated with realizable physical 
situations, and satisfy minimum uncertainty relations. However, those relations only refer 
to the product of uncertainties for pairs of conjugate variables, and not individual uncer- 
tainties, so given some set of {{Oa,^ai5a)}a=i,...,A, not all coherent states may satisfy our 
definition of semiclassical states. To look at this question, let us summarize the freedom 
available in defining coherent states for a theory. 

Consider a theory with a linear phase space F = R^^ = {(gi,Pj)}j=i,...Ar, where the phase 
space coordinates are dimensionless. Then a coherent state centered at {qi,Pi) G F in the 
p-representation is (the square root of) a Gaussian wave function jl5l |. 

ipp,qAp) = njli -^i {-{Pi - Pif/^(^1 - i qi {Pi -Pi)}, (3.4) 

which depends on the N parameters (jj, the widths of the Gaussian in the various Pi direc- 
tions; J\fi = (A/27r(Ti)~^/^ is a normalization factor. The values of the cxj are the freedom we 
have (once we fix the classical phase space point {qi,Pi)) to ensure that the state is semiclas- 
sical. Since the Fourier transform of a Gaussian of width ai is another Gaussian, of width 
l/2(Ti, we expect that for each choice of observables and tolerances {{Oa,€a,Sa)} there will 
be a finite range of values of the for which the state 4'p,g,a is semiclassical, provided the 
tolerances are large enough to be consistent with the uncertainty relations. 
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In the familiar case of a simple harmonic oscillator, where the Hamiltonian contains a 



pair of parameters {mi,Ui) for each degree of freedom, one normally takes ai = ^miUJih/2, 
and ends up with the coherent states commonly used in that case, while other choices for 
the (Tj lead to squeezed states. Here, we don't have such dimensionful parameters available, 
and it might seem that this fact will make the process of defining coherent states somewhat 
difficult. In reality, this is not the case; we don't have the same natural choice for the Uj, 
but we have instead the freedom of choosing values that will make our states satisfy our 
semiclassicality conditions. 

The situation for Bianchi I models is then as follows. Our observables Oa are the p± and 
v±, and we have a single constraint H; since the two quantities in each of the two observable 
pairs have a similar physical interpretation, we will assign the same tolerances to them. Thus, 
given a classical state {pi, Pi) G T and a set of tolerance parameters {ep,ey,Sp,Sy,eH,SH), a 
quantum state ijj will be called semiclassical if it satisfies the conditions 

\{p±)i> - P±\ < (^P , - v±\ < , (Api)^ < (5p , {Av±)^<6y, (3.5) 

and, if ip is not a solution of the quantum constraint. 



Notice that we could have chosen the observables po and vq as part of our set, and assigned 
tolerances to them as well. In this model, our results would not have been very different, 
although in general one cannot even assume that observables which are functions of other 
operators that have small quantum fluctuations have small fluctuations themselves. 

In the kinematical Hilbert space, based on a linear conflguration space, we can use coher- 
ent states of the form ()3.4|) . Then, given a set of tolerances we will have more inequalities 
than free parameters (the or a set of equivalent ones); although it is reasonable to expect 
that solutions do exist for sufficiently large {ea,Sa), one of our goals is to check that this 
is true for "small" tolerances as well. In the physical Hilbert space these states cannot 
be used directly, but we will follow the group averaging procedure to set up appropriately 
modified coherent states peaked around points of the reduced phase space, check if they are 
semiclassical as well, and whether the original kinematical states can be considered good 
approximations. In the concluding section, we will mention models for which the configura- 
tion space is not linear, in which case one needs to replace the coherent states defined above 
by a different choice. The specific form of the states used will depend on the characteristics 
of the model, and we will discuss some options later. 

IV. KINEMATICAL COHERENT STATES 

The first step in the quantization procedure is to find a representation of the CCR on the 
relevant Hilbert space. In our case, given the form of the Harniltonian, the most convenient 
choice for the kinematical Hilbert space is to set Tikin = L^(C,dpj). In this representation, 
the pj operators act by multiplication and the f3i act as f3i = id/dpi] this is the standard 
Schrodinger representation of the CCR. As for the additional observables Vi, since the mea- 
sure /i on Tikin is trivial, we simply get div^ v± = div^ vq = 0, and 




{H)^\<eH, {AH)^<6h. 



(3.6) 





(4.1) 
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where di := d/dpi. 

Given a point P = {Pi, pi) on the constraint surface T in the classical phase space T, 
we construct a coherent state centered at P in the above representation with the usual 
prescription (|3.4|) . as 

ipkmiPi) = Yl.^^^^Afiexp {-(pi - pif/2Ti -i{pi- pi)f3i^ , (4.2) 

where the free parameters Tj measure the width of the coherent state in momentum space 
(in the natural Planck units), analogously to the cij in (|3.4p . and the normalization factors 
are Mi = (ttTj)"^''^. In this state, we recover the usual properties of coherent states, namely 
that the expectation values of the canonical variables are 

(A)kin = A , {Pi)km = Pi , (4.3) 

and their fluctuations have the following expressions 

(AA)L = 1/27:. , (Ap,)L = n/2 , (4.4) 

so that (A/9j)(Apj) = | for each i (no summation over i). In practice, given the similar 
interpretation of /3+ and we will often set r_|_ = r_, for simplicity. For the moment, we 
shall leave the three unspecified, but we will assume that the state is narrowly peaked in 
the sense that the Ti are small, Ti/p^ -C 1 for all i. 

Although we are not imposing the quantum constraint on this state, since we will compare 
it with a corresponding physical coherent state, the quantities we need to calculate are the 
expectation values and fluctuations of the physical operators we have identified, in the state 
■^kin- We already have those of p±; for the remaining ones (among which we include vo, 
although strictly speaking it is not necessary), the expectation values are 

(^^±)kin = J d^P iPtiniPi) i (Po 9± +p±do) ^km(Pi) = Po/3± + P±/3o , (4.5) 
(5o)kin = J d^P ^kin(Pi) i iP+ d- - P- d+) ?Akin(p*) = - , (4.6) 

and for the fluctuations (AOa)^;^ = (O^)kin — (Oa)^;^ we obtain 

(A^±)L = i + i + ^) + I (f + I) + ^ il3lr± + Piro) , (4.7) 



(A^o)L = -I + 1 + ^) + I + - ) +\CP'-n + Plr-) . (4.8) 



In addition, for the Hamiltonian constraint operator we get 

(^)ki„ = H-^o + r+ + T_) (4.9) 
(A^)kin = k iPl^o + PW + v'-r-) + I {rl + ^ + ) , (4.10) 

where we have used the classical constraint — Po+P++P- = 0. Thus, as expected on general 
grounds 0], the expectation values of all relevant operators agree with the classical values, 
except for the (quadratic) Hamiltonian. In addition, for the kinematical coherent state to 
be a good semiclassical state according to our choice of tolerances, the above fluctuations 
need to satisfy the conditions ()3.5|) and ()3.6|) . 
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Allowed values for Kinematical States 




0.2 0.4 0.6 0.8 1 



FIG. 1: The region in the parameter space of (tq, t) where the inequahties for a semiclassical state 
are satisfied is shown together with several of the curves that bound the region. 



Specifically, the only nontrivial condition arising from the expectation values of observ- 
ables is the one from the constraint in ()3.6|) . 

|T+ + r_-ro| <4ej^, (4.11) 

and the conditions on the uncertainties in the p±, the -O-i- and H give, respectively 

r±<2 5l (4.12) 

i + i + ^) + I (f + I) + I (/5oV. + /^iro) < 51 (4.13) 

\ {pU + PW + pV_) + I {tI + rl + T^_)<5l (4.14) 

These inequalities could certainly always be satisfied if all tolerances were considered 
as free parameters, for example by choosing large enough values of the 5a, but what we 
are interested in is finding, for fixed {{ea,5a)}, a range of values for the which satisfy 
the inequalities. We will show explicitly how to do this in the case r_|_ = r_ =: r, which 
simplifies both the calculations and the visualization of the results. 

For example, Eq (j4.14j) becomes now a simple quadratic relation between r and tq, 

tI + ApIto + {ApIt + 2t^-S51)<Q, (4.15) 

which can also be written in a form showing explicitly that (r, tq) must be inside an ellipse 
centered at {—pi, — 2pq). The full set of inequalities is 

0<r<25j (4.16) 
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2r - Aen < tq < 2r + 4e 



1^0 



2(r 



H 



iP^T + PIto) Tot + (1 - 2 51) t^t + \ (r^ + r^) + (p^r + pIt^) < 



(4.17) 
(4.18) 
(4.19) 



where the last inequahty is equivalent to ()4.13|) for no n- vanishing r and tq. 

Fig. [T] illustrates roughly the range of allowed (r, tq) values, for the choice of parameter 
values pQ = 1, p± = l/v/2, (3q = 1, /3± = 1/^2, 5p = 0.2, 6y = 10, en = 0.1, Sh = 0.6. 
Strictly speaking, many values in the "allowed region" shown do not satisfy the conditions 
under which our small-r/pg approximation is valid, but the plot gives a useful qualitative 
description of the effect of each inequality in the tq-t plane. 

The following argument tells us that the first three inequalities can in practice be replaced 
by one of them. Notice that, given that we are assuming pQ to be large, and the relationship 
between H and the p^'s, it is reasonable to set for simplicity 



so 



51 < en, 



(4.20) 



Consider first the inequality ()4.17|) . It defines a strip of allowed (t, tq) values around the 
To = 2 r line, and this strip intersects the r and tq axes, respectively, at 



2 en 



TQ = AeH 



(4.21) 



With the choice of parameters we just made, 2eH > 2 5^, so the first of the two inequalities 
in ()4.17p becomes redundant, provided that ()4.16|) is satisfied. However, if we now look at 
fj4.18|) . it is easy to see that the ellipse intersects the axes at, respectively. 



T=Jpt + ASl 



Pi 



^0 



Apt + 85j,-2pt,<AeH 



(4.22) 



where we have used the fact that pQ is large in the inequalities. Thus, kinematical coherent 
states are semiclassical if (r, tq) lie inside the portion of the ellipse ()4.18p that defines a 
neighborhood of the origin in the first quadrant. 



< To < 2pl 



1 + 



25| 
Pt 



0<t<pI 



+ 4:PlTo-8 5jj 
2PI 



(4.23) 



and satisfy the single pair of additional inequalities ()4.19p . 

We now need to find out whether r and tq can be chosen to meet these requirements. 
Clearly, not all points in the region ()4.23|) satisfy ()4.19p . This can be seen setting tq = 0, 
in which case we get + P±t < 0, which has no positive solutions for r; in fact, if we 
consider ()4.19|1 as a second-order inequality for r for a given (positive) value of tq. 



CpIto + \)r^ + W±rl + (1 - 2 5l) Tq +pI]t+ {\ + plr^) < 



(4.24) 



we see that there are positive solutions for r only if the coefficient of the linear term is a 
sufficiently large negative number. We can also see directly from ()4.13p that we must have 
51 ^ 1. For any positive tq, as 6^ becomes very large, the range of values of r which satisfy 
(|4.24p becomes the full positive half-line (0, oo), so all we need to do is pick a large enough 
6y for all conditions to be satisfied in a region of the t-tq plane. The fact that 5„ must be 
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large for suitable states to exist is an example of the uncertainty relations putting bounds 
on products of tolerances, and a consequence of the fact that we are requiring 5p to be small. 

We finish this section by giving a set of parameter values (different from the ones in Fig.^ 
which satisfy the relationships mentioned above, and allow us to choose smaller values for 
r and tq, for use in later parts of the paper. If for simplicity we set p+ = p_ and /5+ = /3_ 
(the two inequalities (j4.24|) or their equivalents become a single one), it is natural to make 
the choices 

p± = l, Po = V2, ^± = 1, ;3o = l (4.25) 

for the classical phase space point (notice that all inequalities are invariant under an appro- 
priate rescaling of the variables involved; this can be used for example to assign any desired 
value to po; which we think of as an arbitrary 'global scale'), and 

5p = 0.15 , 5^ = 12, eH = 5H = 0.25 (4.26) 

for the relevant tolerances. Then, a pair of width parameters that solves the inequalities is 

r = To = 0.025 . (4.27) 



V. GROUP-AVERAGED COHERENT STATES 

In the Dirac approach to the quantization of first-class constrained systems, a physical state 
"^phy is one which is annihilated by the Hamiltonian constraint, H ■ "^phy = 0. It is a general 
feature of constrained systems that physical states will not belong to the kinematical Hilbert 
space Tikin, and the Refined Algebraic Quantization program has been developed to deal 
with such instances [16]. The general strategy is the following. One first identifies a dense 
subspace $ of the Hilbert space Hkin, where the Hamiltonian is well defined (that is, $ is 
contained in the domain of H). The physical states will then belong to the dual $*, which 
satisfies the inclusion relation $ C Tikin C $*. The choice of the dense subset $ is a delicate 
matter, and normally has some physical implications [3], but for our purposes we are only 
interested in looking at coherent states which have some "nice" properties that make them 
well behaved, and those states belong to most choices of the space The procedure we will 
follow for implementing the Dirac requirement that physical quantum states be annihilated 
by the constraint operators is the so-called group averaging, a particularly popular approach 
within the general RAQ program. 

The intuitive idea is to start with a state and average it over the group generated 

by the constraint. If all goes well, one will end up with an element of $* satisfying the 
constraint. Our strategy will be very naive. Starting with the kinematical coherent state 
(j4.2|) . group averaging gives us the physical coherent state ■i/'phy = A^~^ / dA e"''^-'^ "ipkin, where 
iV is a normalization factor. Since in pj-space the constraint is a multiplication operator, the 
integration is straightforward and gives a delta function 6{pl — — p^), so we can identify 
'0phy(p±) with tpkm{{p+ + P-Y^'^ ,P±), up to a new normalization factor to be determined,^ 

i'phy{p±) = A/'e"t(P++P-)'^'~P"]''/^^°"(P+"^+)'/^^+"(^'~"^-)''/^^- X 

^ g-i[{p++P?.)^/^-Po]/3o-i(p+-p+)/3+~i{P--P-)/3- l\ 



The group average procedure, in a strict sense, will yield also a contribution from the past null cone, but 
since we are considering Gaussian wave functions that are almost zero there, the contribution from that 
null cone is negligible. 
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Of course, strictly speaking what one is doing is defining a function on the future null cone 
{Pq — P+ — P- = 0, Po > 0), with coordinates p±. The inner product induced by the group 
averaging procedure is then of the form 

^'^''^^ " / 2 {pl'+^p^y/^ ^phy(P±) (l^phy{p±) . (5.2) 

Again, we want to calculate the expectation values and fluctuations of the physical observ- 
ables p± and Vi, and compare them to the values computed with the kinematical coher- 
ent states. The strategy here will be to represent the Dirac observables in the physical 
Hilbert space making use of the measure ()5.2|1 induced by the group average and the gen- 
eral form of the operators given by ()2.3|1 . The momentum observables act as multiplication 
operators; to find the action of the configuration-like ones we use again ()2.3j) . where now 
IJ,{p±) = \ {p\ +p'i)~^/'^ is the measure appearing in ()5.2|) . This gives, on physical states. 



v± ip{p±) = i^Jpl+P- d±il){p±) , V(p±) = i {p+d^ - p^d+) iIj{p±) . (5.3) 

The first step, before calculating expectation values and fluctuations, is to calculate A/". 
However, because of the presence of the expression (p^ +p^)^/^ in the exponent of ipphy, as 
well as in the measure, the relevant integrals cannot be calculated exactly, and we need to 
use an approximation. 

If the state -^phy is narrowly peaked in the p± directions, we can make use of the assump- 
tion r -C Pq in the approximation, and keep only the first few terms in the expansion of 
every result in powers of t/pq. To see how to implement this idea, let us parametrize the 
deviation of p± from the classical values p± by a new pair of variables, 

X := — cos ^ + — sin ^ — 1 , y := — — sin ^ + — cos 6 , (5.4) 

Po Po Po Po 

where 6 is the angle between the vector (p+,p_) and the positive p+-axis, and the origin 
in the (x, y) coordinates corresponds to the point where the Gaussian is peaked. We can, 
equivalently solve 

p+ = {x + l)p+ - yp^ , p- = {x + l)p^ + yp+ . (5.5) 

Intuitively, x and y are the deviations from p± along the radial direction in the p± plane 
and along the + = pi circle, respectively, both expressed in units of po- Using these 
definitions we get, e.g., that p\+p'^ = pi [(1 + x)^ + y"^], and in terms of the variables x and 
y the physical state can be written as 

ijp^yix,y) = Are-^[(i+^)'+J^']'''-i}'^^o/2ro-(:.^+y^)p1/2r ^ 

in the remainder of this paper, we will assume for simplicity that r_|_ = r_ =: r. Notice also 
that we will always assume that the parameters satisfy the classical constraint pi = p\+p^_. 

Our strategy for obtaining expressions for A/", as well as all expectation values and fluc- 
tuations, will be the following. All calculations can be reduced to integrals of the form 

'^^'^^ i'{x,yy 0^{x,y) , (5.7) 



[(l + x)2+|/2]l/2 
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where O may be a combination of multiplication and derivative operators. Since ijj*Oijj 
always contains the Gaussian factor exp{ — (x^ + y'^)pQ/T} which is peaked at small values 
of X and y, we will expand all remaining factors in the integrand, including the measure, 
in powers of x and y, and calculate the integral term by term in the expansion. Each 
non- vanishing term is a Gaussian integral of the form 



dxa;2"e-"'''oA 



dyy 



2m 



^ 7r(2n-l)!!(2m-l)!! [ t_ 

2^n-\-m 



n+m+1 



Po, 



(5.8) 



with a coefficient that depends on the parameters in ip and, most importantly, contains one 
inverse power of r for each derivative in O; if we want to keep terms up to a certain order in 
t/pq when adding the integrals of the type (j5.8|) . we need to take this into account. With the 
exception of {AH)1^^ and {Avi)^^^, the expectation values and fluctuations of our operators 
in the kinematical states contain a single power of r; therefore, in order to compare those 
results with the corresponding ones for physical states, for the latter in practice we need to 
keep terms in our expansions up to the two leading orders in r/pl- 

The normalization factor J\f is the simplest expression to calculate. Since 



dx dy 



[(l+x)2+y2]l/2 



|V^phy(x,y)P = ^ 



r 



{5S) 



the result is 



2po 
vrr 



1 + 



'pI 



V 



T 
"i^2 



T 



(5.10) 



To -zj 2p^o_ 

With Af, we can now calculate other expectation values. It is useful to first calculate the 
auxiliary quantities 



2fo 



T 



- 1 



'pj_; 

.^0 



4p; 



>^4 



0(r3 



{y) = o 



T 



with which we find for the p±, 



{P±) 



phy 



dp+dp_ 



2(p2_+p2)l/2 

P± {{x) + 1) =p± 



p± 



phy I 



.^0 



(5.12) 



and for the fiuctuations of the same operators. 



(P±)phy 



pl{l + 2{x) + {x^))+pl{f) 



-2 _^P% 
P± + — 



pI 



2fo 



P% 
Api 

T 



T 



(^P±)5hy = (P|)phy - (P±)Jhy = ^ + 



'1 

2 To, 



2pl 



+ 0(r=^ 



(5.13) 
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As for the Hamiltonian constraint, the exact physical state ?/'phy of course gives {H)phy = 
and (A_f/')pj^y = 0, by definition. 

The Vi operators contain derivatives with respect to p±] to calculate their expectation 
values it is useful to recast those in terms of {x,y), and we find, from ()5.3|) and ()5.5p . 

^ \Po Po J 

vo^p{x,y) = i[ - yd^ + {1 + x)dy]ilj . (5.14) 

Then for the expectation values we find 

= {PoP± + P±Po) - {Poh + 2p±^o) ^ + 0{t^) , 

(^o)phy = '-j ^(^ ^ %^ly2Yl2 ^phy ( ' V 9. + {I + x) .9,) 

= {pj. -p.p^) (l - ^ + 0{r')^ ; (5.15) 
similarly, the fluctuations of these operators are, respectively, 

^^"^W=27 + 4(^-^ + — 

i^^o)l^y = §-\ + 0{r). (5.16) 

We now have all the results we need to write down the conditions that the physical 
states be semiclassical, but we can also bypass this step by a direct comparison between the 
expectation values for the Oa in the two sets of states, 

(p±)phy-p± = -^ + 0(r2) (5^17) 

{v±Uy -v± = -iPoPl + 2p±^o') ^ + 0{t^) , (5.18) 

and between the corresponding fluctuations, 

(Ap.&,-(Ap.)L^-(f + + (5.19) 

(A«j;^, - (ACJL = J fi + 11 + L\ - - j3lr„ + 0(t) . (5.20) 

4 \ T To / 4 Po 

The first pair of relationships tell us that the physical coherent states have expectation values 
for the Oa that are close (to order r) to the classical values; that fact, together with the 
observation that the right-hand sides of the last two equations are negative (provided that 
To < 2pq), shows that the physical states are also semiclassical, with respect to the same 
tolerances used for the kinematical state and values for and e„ that can be chosen to be 
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small of order r. In particular, if we choose as reference classical state the one specified by 
fl4.25|) . and our criteria for a semiclassical state include the choices ()4.26p plus, for example. 



ep = 0.01 , = 0.02 , (5.21) 

then the state ipphy constructed using the width parameters given in ()4.27|) is semiclassical. 

Thus, with this choice of widths, the kinematical state is compatible with the physical 
one in the less restrictive of the two possible definitions mentioned in Sec IIIII In the more 
restrictive sense, however, it is in general not a good approximation, because the right-hand 
side of (|5.19|) is not necessarily small compared to (Ap-i-)pj^y, nor that of (|5.20|) compared to 
(A'O-i-)pj^y. For the former to happen, the widths would have to satisfy r <C tq (or p± would 
have to be very small, but this cannot be true for both signs), while for the latter to happen, 
they would have to satisfy (either r ^ tq or) r ^ tq ^ Po- 

VI. TIME EVOLUTION 

In the previous sections we have considered two kinds of coherent states, kinematical and 
"group averaged" ones. The first kind are coherent states defined on the full phase space 
and peaked around a point of the constraint surface. They do not satisfy the quantum 
constraint. The second kind of states are constructed out of the first kind by averaging over 
the orbit generated by the constraint. Even though the second kind of states are in a sense 
dynamical since one can think of the constraint (for totally constrained systems, such as this 
one) as generating time evolution, in a strict sense they are 'frozen'. Once the quantum state 
belongs to the physical Hilbert space, the notion of 'gauge' time evolution is lost. There is 
no notion of time evolution, no dynamics. Within this perspective, the observables that we 
have considered, namely Dirac observables, are constant along the orbits of the constraint 
and are therefore constants of the motion according to the same dynamical interpretation. 
This is the well known "frozen dynamics formalism" for time-reparametrization invariant 
constrained systems. 

Is there a way of recovering dynamics from such a formalism? Again the answer is in the 
affirmative and several solutions have been known for a while jiol . [l7| . We shall consider 
one of those possibilities. For this we shall proceed in two steps. In the first part of the 
remainder of this section we recast the quantum constraint equation in such a way that 
it resembles a Schrodinger equation. In the second part we consider observables (different 
from the ones already considered) that capture the intuitive notion of time evolution and 
that will allow us to ask dynamically meaningful questions such as: given a coherent state 
with a small spread in the anisotropics at, say, the Plank time, we would like to know how 
much time later (as measured by a proper measure of the total volume) the state will spread. 
Based on the previous results obtained so far, where we found a range of parameters that 
produce acceptable coherent states, we shall try to give answers to these type of questions. 

A. Schrodinger Equation 

The idea here is to consider states at the kinematical level that are in the p representation 
for the ± sector but are functions of the spatial volume and will therefore be diagonal on /3o. 
That is, we will consider functions of the form ^I/(/5o,p±). We want to consider the quantum 
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constraint H = — | (Pq + + p^) =0, rewritten as 



H = -\ {Po + ^Pl + )(po - ^Pl +Pl). (6.1) 
and look for solutions to the quantum constraint that are of the form 



(po-Vn + r-)-^(/5o,p±) = o. (6.2) 

The choice of relative sign in (|6.2j) is equivalent to the choice we made earlier of the future 
(as opposed to past) light cone in p-space, in the context of group-averaged coherent states. 
In a sense we are following the original "already parametrized" program of ADM. 

Note that this strategy is equivalent to considering the reduced phase space, described by 
true Dirac observables, and defining wave functions of them that have a "time dependence" 
as they are also functions of [3q. Among these states, we can look for coherent states peaked 
around 'physical' phase space points. We start therefore with the reduced phase space 
Fr for which we will use the coordinates = {(p/,/?/), / = ±}, and the Hilbert space 
Hr = L^(1R^, dp+dp_). The sector of the solutions to the Hamiltonian constraint now takes 
the form of a Schrodinger equation, 

. d 



^(/?o,P±) = Vp+ ^iPo,P±) . (6.3) 



That is, any (normalizable) function ip{p±) gives rise to a solution to the quantum constraint 
via the "time evolution equation" ()6.3|) . whose solution can easily be written as 



^(/5o,P±) = e-'V^Po^{p^) . (6.4) 

The strategy now is to start with a coherent state "intrinsic" to the p± plane (as a 
coordinatization of the reduced configuration space) and then construct the corresponding 
physical state. This choice is motivated by the fact that the quantities P± have a clear space- 
time interpretation as the anisotropics of the cosmological model. Let us then consider, as 
initial state, 

V'o(p±) = n exp[-(pj - Pif/2Ti -i{pi- Pi)(3i] , (6.5) 
i=± 

which is a coherent state peaked around the point (p±, /5±) on the reduced phase space, and 
with spread dictated by the parameters Xj. We get then 

^o(/5o,P±) = e-'^/^-^'MP±) 



n M exp [-(p, - p,)V2r, - i (p, - pM e-^V^*. (6.6) 



i=± 



Let us note that the general solution should be a function of (/?o ~ ^0)5 for ''^o the "initial 
value" of the parameter /Jq. Now, when (3q = 0, the classical volume element is given by 
y/det q = 1 in Planck units. Thus, the "Planck epoch" is given by the values of /?o = (recall 
that the "singularity" is at (3q — * —00). We shall then, in the remainder of this section, take 
to = 0, and consider the "initial wave function" at to = as being defined at the Planck 
epoch. 
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B. Time Evolution of Observables 



The next step is to consider observables. On the one hand we know from the general 
theory of constrained systems that physical observables are those operators that preserve the 
space of physical states. Properly represented Dirac observables have that property and that 
is why we have considered them in previous sections. On the other hand, we are interested in 
disentangling dynamical information from the particular (Schrodinger-like) representation 
of this part. This means that we need observables that capture the notion of time evolution. 
Such observables exist and are sometimes referred to as evolving constants of the motion 



17[. The physical states considered in this part are solutions to the equation ()6.3|) so we 
have to define observables that preserve the space of such solutions. If we were to consider 
the operator /3o (that naturally acts as a multiplication operator), we would be thrown out of 
the physical Hilbert space. It is not a physical observable. What we need to do is to consider 
the following situation: let us assume that the phase space function (3q is an (internal) time 
parameter, and that we can consider a new kind of observables 'defined at time' /3o = t. 
That is, we fix a /5o = t slice in the configuration space {p±, Pq) and consider the restriction 
of the wave function to that slice, "^{Pq = t,p±). We can then act with an operator (that 
might be '/5o-dependent') and consider the time evolution of the resulting wave function, 
thus obtaining again a physical state. This is nothing but the ordinary prescription for 
defining Heisenberg observables in ordinary quantum mechanics. 

Consider an observable that depends explicitly on /Sq and does not commute with the 
constraint (for instance the volume form V{l3o) = e^^°). The strategy is to define a one- 
parameter family of observables V{t) = V{i3q = t) as Heisenberg observables and then 
consider the corresponding operators. Note that V{ti) ^ V{t2) if ti ^ t2, as physical 
operators on the physical Hilbert space. The observables that we will be interested in are 
of two kinds. On the one hand, we shall consider the 'configuration observables' p± that 
are constants of the motion since they commute with H, and on the other hand those 
operators that measure the anisotropics, namely the operators /3±. These are not explicitly 
'time dependent', but we expect both expectation values and dispersions to be changing 
quantities.^ 

Let us then compute the expectation values and quantum fluctuations for the operators 
/3± = id/dp±. The expectation value will take the form 



0^)t = 



(^p± e" 



P+ 



p2_ 



t + (3, 



(6.7) 



Note that there will be a shift of the 'peak' of the wave function, since the expectation 
value will have a linear t dependence. This is to be expected since classically we know that 
the system evolves in a way that the f3± depend linearly on time. In order to compute the 
coefficient we need to make some approximations as in previous sections. Making use of the 
same assumptions, we arrive for the integral to the following expression. 



0+)t = P+ + tM^pl / dxdy e-pl^-'+y'y 



[I + x)p^ 



yp~ 



Po[(l + a;)2 + y2]i/2 



(6.8) 



^ Let us qualify the statement: We expect that the one-parameter family of numbers {(3±)t, corresponding 
to the one-parameter family of operators /3± will depend on the value of t. 
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This can be approximated, as a power series expansion in r, by the expression 



Po 



1 ^ 3 



(6.9) 



Note that at zeroth order in r the speed of the peak coincides precisely with the classical 
value. Let us now compute the fluctuations of the operators and consider the expectation 
value to second order in r as, 



0l)t 



2t ^ Po 



1 r 
4M 



3 r2 



t + 



pi 
pI 



It 33 r2 
2fo^^¥o 



(6.10) 



From here we can then compute the fluctuations of the anisotropy observables. We have 
then, 



0l)t - 0±)l 



1 143 pi 2 2 
2t 16 



(6.11) 



where we keep terms that are at most quadratic in r. We see that, as in the case of a 
non-relativistic free particle, the spread of the Gaussian grows quadratically with 'time' t. 
The next step is to estimate the value tmax of the maximum allowed value for t before the 
wave packet spreads considerably. For that we use the values of the parameters p±, po and 
T chosen in (j4.25|) and (|4.27j) . The reason for this is that we expect that the conditions 
imposed on the expectation values and fluctuations of the relevant operators yield realistic 
values for the parameters defining the coherent states; even though in this part we start 
with intrinsic coherent states at 'Planck time', the resulting wave function is closely related 
to the physical coherent states of Sec. El and it seems reasonable to employ the parameters 
used there. If we define the time for the wave packet to spread significantly as the value 
for tinax such that {Ap±)f^^^ is about twice as large as {Ap±)^^ ^ 20, with the values of the 
parameters given in Sec IIVI it is easy to see that t^ax ~ 170. That is, the time for which 
the packet spreads to about twice its size is within 200 times the fundamental time unit, 
the Planck time. With a different choice of parameters that yield smaller values of r, one 
might have slightly longer spread times, but within the same order of magnitude. Another 
possibility would be to take the values used in Fig. [T]and choose r = .03, which is near the 
high end of the range of allowed values. In this case the time tmax ~ 65, again of the same 
order of magnitude as with the previous choice. 

One might ask about the possibility of having semiclassical states that approximate the 
/3's more accurately, which means smaller initial fluctuations (of the order of 1/t). This 
would require large values of r, for which the expansion in powers of r would not be justified. 
In that case, one would need to expand expectation values and fluctuations of observables in 
inverse powers of r, and use those in the set of inequalities that express the semi-classicality 
conditions. With this choice, one might be able to construct states that take longer to spread. 
We shall not explore this possibility further. Let us end this section with a remark. Even 
though we have found a 'time' (as measured in terms of the volume element) for which the 
initial wave-packet spreads, we would not like to suggest that this has a direct significance 
for the spread of the wave-function in realistic cosmological models, since we have neglected 
any matter and the effect that this may have in possible decoherence effects on the wave 
function. 
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VII. CONCLUSION AND OUTLOOK 



In this paper we have only considered the simplest type of anisotropic cosmological models, 
the vacuum ones of Bianchi type I. For our purposes, the simplicity of these models lies 
in the fact that their phase space is a vector space, their only constraint is quadratic, and 
it depends on the p variables only; in this sense, these models are similar to the general 
class of models treated in Ref 0], in which kinematical and physical coherent states are also 
discussed. However, on the one hand our Hamiltonian and the set of physical observables we 
chose are not of the type used there, and on the other hand the fact that the Hamiltonian 
itself is the constraint raises the usual issue of time evolution in generally covariant systems, 
and led us to the construction of a third type of state that we called "intrinsic", with 
which physical predictions on the evolution of physical observables can be attempted. The 
relevance of our findings in this simple model for the ultimate goal of constructing realistic 
semi-classical states in full quantum gravity still needs to be explored. However, we hope 
that these first steps might shed some light on the program. 

To explore further the relationships between kinematical and physical coherent states as 
candidate semiclassical states for gravity models, a natural extension of this work would 
be to consider other Bianchi type models. As soon as one starts doing that, one is faced 
with the fact that the potential f/(/5*) does not vanish, which implies that the quantum 
Hamiltonian is no longer a multiplication operator. Fortunately, in some cases this difficulty 
can be circumvented. As is well known, the dynamics of a system on a configuration space 
with a given metric subject to a potential U, can be mapped into that of a system with a 
new metric (the Jacobi metric) and vanishing potential. In our case, the transformation of a 
Hamiltonian like (j2.2|) into an equivalent one describing a free particle follows the standard 
procedure, but what one usually ends up with is a curved Jacobi metric replacing rjij. What 
is not so obvious is the fact that in the case of many homogeneous cosmologies, known as 
diagonal, intrinsically multiply transitive (DIMT) models, the new metric is still fiat jlol |. 
These models include the following cases: Bianchi types I and II; sub-families of Bianchi 
types III, VIII and IX defined by /3~ = (the Taub model); Kantowski-Sachs spacetimes; 
and Bianchi type V models with f3^ = 0. For Bianchi type I and II, the minisuperspaces 
are 3- dimensional, parametrized by f]~^ and (3~ in the Misner scheme; for the remaining 



DIMT models they are 2-dimensional, parametrized by an appropriate subset of the /3* [11 
Note that even though the type HI and V models are class B, which means that in general 
the Hamiltonian description is more subtle fI2], these particular cases do not suffer from the 
problems described in Refs [llll^. 

Thus, DIMT models are natural candidates for a first extension of our work. When 
one carries out the transformation to a new set of phase space variables (/9*,Pi) in which 
the potential term in ()2.2|) vanishes and the metric again has the form rjij, the difference 
between the various DIMT models shows up in the fact that the allowed ranges of the pi 
are not the whole real line 0|. Therefore, those variables cannot be regarded as momenta 
in the regular sense of elements^ of the cotangent space, and in the quantum theory the 
corresponding operators in the /5j-representation cannot be the usual derivative operators. 
(Note: this feature, which also occurs in full quantum gravity, has been addressed in that 
context by Klauder 1^ and LollJS.) We have then one more reason to set up the quantum 
theory in the p-representation [lO| . This time, however, the "configuration space" for these 
models is a proper subset of 2+1 Minkowski space, with restrictions on the values of the pi 
depending on the model under consideration, and the constraint surface is the intersection 
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of the (future) light cone of the origin with this subset. 

With this setup, one can define suitable generalizations of coherent states in the config- 
uration spaces of the various DIMT models, which now are non-trivial spaces (for example, 
a half-line rather than the whole real line), and analyze their properties as candidate semi- 
classical states. One possible strategy is to construct such states as suggested by Klauder's 
work ji^l, and search for a proper subset of states that are consistent with our physical re- 
quirements, but we shall leave that investigation for a future publication. As a further step, 
one can envision extending the work to models which cannot be formulated as free particles 
on portions of Minkowski space, which will require new techniques to handle the quantum 
constraint (it would be interesting, for example, to study the connection with Kiefer's use of 
the "principle of constructive interference" to build wave packets in minisuperspace |2ll an d 
other models for the emergence of classical Friedmann- Robertson- Walker spacetimes j22|). 
and then possibly some non-homogeneous, midisuperspace models. In terms of our original 
goal of understanding the relationship between the kinematical and physical semiclassical 
states for quantum gravity, the hope is that at some point one may see a pattern that will 
allow us to make statements of a more general nature. 
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